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Blandford-Znajek mechanism, by which the rotational energy of a black hole is extracted through 
electromagnetic fields, is one of the promising candidates as an essential process of the central engine 
of active compact objects such as Gamma-Ray Bursts. The only known analytical solution of this 
mechanism is the perturbative monopole solution for Kerr parameter a up to the second order terms. 
In order to apply Blandford-Znajek mechanism to rapidly rotating black holes, we try to obtain the 
perturbation solution up to the fourth order. As a result, we find that the fourth order terms of the 
vector potential diverge at infinity, which implies that the perturbation approach breaks down at 
large distance from the black hole. Although there are some uncertainties about the solution due to 
the unknown boundary condition at infinity for the fourth order terms, we can derive the evaluation 
of the total energy fiux extracted from the black hole up to fourth order of a without any ambiguity. 
Further more, from the comparison between the numerical solution that is valid for < a < 1 and 
the fourth order solution, we find that the fourth order solution reproduces the numerical result 
better than the second order solution. At the same time, since the fourth order solution does not 
match well with numerical result at large a, we conclude that more higher order terms are required 
to reproduce the numerical result. 

PACS numbers: 04.70.-s,95.30.Qd,95.30.Sf,97.60.Lf 



I. INTRODUCTION 

Gamma ray bursts(GRBs) are the most powerful ex- 
plosions in the Universe [J, Q. GRBs release 10^^ — 
lO^^erg in a few seconds, which corresponds to a rest 
mass energy of the sun. The mechanism of the central 
engine of GRBs is not still understood well. However, 
from the fact that some GRBs are associated with su- 
pernova explosions 0, Q , we can guess that there might 
be a black hole at the center of a GRB and its energy 
source may be the rotation energy of the black hole. In 
order to extract the rotational energy of the black hole, 
the Blandford-Znajek mechanism was proposed [EBQ- 
In this mechanism, the rotational energy is extracted in 
a form of the Poynting flux that is realized by the sur- 
rounding electromagnetic fields. The scale of extracted 
energy from the black holes is, in a rough estimate Q, 

E ~ B'^njjM'^ 

where i3 is a magnetic field, M is the mass of the black 
hole, a is the specific angular momentum of black hole 
(Kerr parameter). Compared with the energy scale of 
GRBs, this energy scale seems to be a good one. 

Although Blandford and Znajek (1977) obtained an 
analytical solution for this mechanism by assuming the 
shape of the magnetic field (monopole magnetic field), 
this solution has two points that should be investigated 
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further. First, this solution is the perturbative solution 
for Kerr parameter a. Thus, this solution is valid only 
for a slowly rotating black hole (a ^ M). Second, it is 
very difficult to impose proper boundary conditions. The 
equation we have to solve is the second order differential 
equation, thus we must impose two boundary conditions 
for a solution. In the original paper, assuming force free 
condition in all region, they imposed boundary conditions 
such that at event horizon the magnetic field is finite and 
at infinity the solution should connect with another force 
free solution obtained in Newton gravity Q . The bound- 
ary condition at event horizon is plausible Q, but since 
at large r, there is no guarantee that we can use force 
free condition, we have to be careful for the boundary 
condition at infinity. 

The aim of this paper is to evaluate more accurately 
the energy fiux extracted from a rotating black hole 
by the Blandford-Znajek mechanism, i.e. to derive the 
fourth order terms of Kerr parameters for the monopole 
solution. Although there is two natural boundary con- 
ditions for the second order solution, proper boundary 
conditions are not known for the higher order solutions. 
Thus, when we solve the field equations, we impose a 
boundary condition only at the event horizon while we 
do not impose any boundary condition at large r in this 
study, and we investigate the behavior of the magnetic 
field at large r. In particular, we examine whether we 
can impose the force free condition at infinity in order to 
check the validity of the force free boundary condition. 
Further, we compare the analytical solution we obtain 
in this study with the numerical calculation in order to 
evaluate the contribution of the fourth order terms to the 
extracted energy flux from the black hole. 

The plan of this paper is as follows. In section 2, we 
solve the perturbation equations for Kerr parameter a 
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by imposing a boundary condition at the event horizon 
and estimate the total energy flux extracted from a ro- 
tating black hole. Here, we use Kerr-Schild coordinate 
since this coordinate is regular at the event horizon. Al- 
though the second order solution is shown to be regular 
in a hole region, we show the fourth order solution di- 
verges irrespective of any boundary condition at infinity, 
which means the perturbation approach breaks down at 
large distance from the black hole. Our notation follows 
McKinney and Gammie Q in section 2. In section 3, 
we compare the analytical solution obtained in section 2 
with numerical solution which is valid for < a < M . In 
section 4, we summarize and discuss our results. 



II. METHOD OF CALCULATION 

As stated in section 1, in Blandford-Znajek mech- 
anism, rotational energy of a black hole is extracted 
through electromagnetic fields and this extracted energy 
may be the energy source of GRBs. When the energy is 
extracted through electromagnetic fields, it propagates in 
the form of the Poynting flux vector. Here, we give the 
formula of Poynting flux vector in Kerr space-time for 
arbitrary shape of magnetic fields. As mentioned above, 
we use Kerr-Schild coordinate here for convenience: 



{-i + %)de + {i + %)dr^ + ^de^ 



4ar sin^ 9 
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Aa^ sin^ O) dcf^ + —dtdr 



^ dtd<j), -2a sin^ 9 (l + ^ j drdcj) (2) 



where 



E = -t- cos^ 9 
A = - 2r + 



(3) 
(4) 



and we use the unit in which G — c ~ 
M (black hole mass) = 1. The energy momentum ten- 
sor of electromagnetic fields is written in terms of F^^ — 



T, 



EM 



F.,„FP - -F,BF"f^ 



(5) 



— T'J component of T^, is the r component of Poynting 
flux vector. 

We can expect that matter's contribution to the total 
energy can be ignored near the black hole, so force free 
condition is valid [s*]. Force free condition is represented 
by electromagnetic tensor F^^, as follows: 



the total energy momentum tensor is composed of elec- 
tromagnetic one only: 



total -^matter 



rpfJ,iy 

^EM- 



(7) 



In other word, matter energy is ignorable compared with 
electromagnetic energy. Thus, in the force free region, 
equation of motion of electromagnetic field becomes as 
follows: 



^ t^^EM — 



(8) 



We use T^'' as T^^/ below. 

Generally, the number of independent variables is six; 
electric fields and magnetic fields. However, when we im- 
pose the force free condition, there are four independent 
variables. Further, assuming the the force free condition 
([5]), stationarity, and axisymmetry, there is a relation 
among vector potential terms : 



At,eA^^r = 



(9) 



which makes the number of independent variables three. 
This relation implies that part of the electric field com- 
ponents are represented by magnetic field components. 
Also, we can define the angular velocity JIf of electro- 
magnetic fields from Q: 



At.6 

Ach.f) 



A, 



A 
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All components of the electric fields are represented by 
this angular velocity JIf and magnetic field components 
B"^ {i = r,9,(j)). After all, independent variables are (j) 
component of magnetic field B'^, cj) component of vec- 
tor potential and angular velocity of electromagnetic 
field flF- Here, we define magnetic field components as 
follows: 



B" 



B'P = 




(11) 
(12) 

(13) 



Using these variables, we can evaluate the energy flux 
— T J of electromagnetic field as 

-T"; = -2[B' fVLF {^F - — ) sin^ O-B'' B'^'^fI^ sin^ 9. 

On the event horizon (r = r+ A = 0), this formula 
becomes 



(15) 



Fu,r = 



(6) 



where is the current vector. This approximation 
means that the inertial force by matter is ignored and 



where fin is the angular velocity of the black hole. 
Thus rotational energy can be extracted through elec- 
tromagnetic fields in the force free region if and only if 
< f^F < ^H- 
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What we want to evaluate is the total energy flux of 
the extracted energy from the black hole. Its total energy 
flux E is represented as the integration of the energy flux 
at r = const, surface. 



E 



-27r 



(16) 



From stationarity and axisymmetry of Kerr space time, 
this total energy flux E is independent of r at which the 
integration is operated. 

The formula ()16p which we derive here can be used for 
any shape of magnetic fields. In this section we evaluate 
this total energy flux for the monopole magnetic field. 
Since it is very difficult to solve the equation of motion 
([5]) analytically, we treat Kerr metric as the perturbation 
in Kerr parameter from Schwarzschild metric to obtain 
analytical solutions. 



A. second order terms 

Here, we explain the monopole solution which Bland- 
ford and Znajek obtained by perturbation method for 
Kerr parameter a. At first, the unperturbed solution is 
the monopole solution in Schwarzschild space time. This 
unperturbed solution satisfies 



y(0) = 



(17) 



where (0) means that Schwarzschild metric is used. The 
first order monopole solution is 



A^^^ = -Ccos6', 



= 0, 




(0) 



(18) 
(19) 
(20) 



where C, strength of the magnetic field, is the constant. 
In this solution, since electromagnetic fields do not ro- 
tate, rotational energy of black hole cannot be extracted. 

In order to obtain the monopole solution around a 
slowly rotating black hole, we consider perturbation for 
Kerr parameter a. Taking account of axisymmetry, we 
expand the solution up to the second order as follows: 



= -Ccos6' -f a'' A 



2^(2) 



(21) 
(22) 
(23) 



field equations. Here, we impose boundary conditions. 



A 



(2), 



Ion the event horizon 



finite. 



A 



(2) 



0(-)(r-oo), 



on the event horizon 



finite. 



(24) 
(25) 
(26) 



The conditions ((24|) .(l26 p in Kerr-Schild coordinate corre- 
sponds to the condition in Boyer-Lindquist coordinate 
such that the magnetic field that FIDO (fiducial ob- 
server) feels should be finite [lOj , and (p5|) is required to 
connect the solution to the Michel's solution. In ref [H], 

(2) 

regularity at infinity and separability of are taken as 
the boundary conditions. Although the boundary condi- 
tions are different from ours, the result we obtain below 
is the same with McKinney and Gammie [6]. As we see 
in fourth order calculation, A''^^ is not separable. Thus 
we do not adopt the separability as a boundary condition 
in this study. 

Under these boundary conditions, the solutions are de- 
rived as follows: 



A'-^> = C/(r) cos 61 sin^e, 
4^2 U r / ' 



(1) 



where 

/(0^(-2(Vln(I-^)ln(^))^ 



(27) 
(28) 

(29) 



1 + 3r - 6r^ 
12 
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1 



^^^2^ + 72 + 3. 



r 



and 



Jo t 



(30) 



(31) 



By choosing uj^^") = 1/8, we can use the force free condi- 
tion at large r consistently. We evaluate the total energy 
fiux of the electromagnetic field. 



E = -2tt 



-gT\de 



2^2 



24 



(32) 



Eq. (121D is the result that Blandford and Znajek (1977) 
obtained. 



We also use the expanded Kerr-Schild metric up to the 
second order in a. 

The equation we solve is ^f^T^ = 0. The t, 4> compo- 
nents of this equation represent the conservation equation 
for energy and angular momentum of the electromagnetic 
fields, respectively. The r, 9 components are the trans- 



B. fourth order terms 

Here, we try to obtain the fourth order terms for Kerr 
parameter. Although the method to derive the fourth 
order terms is basically same as the second order case. 
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we have a problem with the boundary condition as men- 
tioned in section 1. We must impose two boundary con- 
ditions. One is the condition at the event horizon. This 
is same as the second order solution. As for another 
boundary condition, we cannot guess anything since we 
do not know where the force free condition breaks down 
and which solution we can connect to. That is why we 
impose a boundary condition only at the event horizon 
and we solve the fourth order equations. We also exam- 
ine whether we can put a force free boundary condition 
at large r using the derived fourth order solution. 

The equation we solve is the fourth order one for Kerr 
parameter a of V^T^ = 0. We expand the variables as 
follows: 

A^ = -Ccose + a'^Af + a^Af , (33) 
B'>' = aB'l'^^'> +a^B'^^^\ (34) 
np = auj^^^ + a^oj^^\ (35) 



By checking the fourth order terms of the t and 6 for the 
equation V^T^ = 0, one can guess that the 9 dependence 
should be as follows: 



Af' = hi (r) cos e + hs (r) cos^ 6" + /15 (r) cos^ 6*, (36) 
S^^'' =5o(r)+ff2(r)cos2^, (37) 
=6o + 62cos26' (38) 



where 60,62 are constants. When we impose boundary 
condition at the event horizon like the second order so- 
lution, we can solve the (f> component of the equation for 
go,g2- The solutions are as follows: 



go{r) = 



92{r) = 



1 / j,4 

——^7 — -288 + 456r + 92r^ + lOQr^ + 36r^ - 11526or3 + ir(67r^ - 521 + 11526o) 

576r^(r — 2) \ 2 

-|-3r^(8 + 22r - ^Qr^ + 6r^ + 12r'^) + 9r^(12 - llr - + 2r^){hiC—-^) In(-) - Li2(-)) 

r 2 r 

1 / 

——5- — -96 + 120r + 20r^ + 125r^ + 24r^ + 12r'' - 38462r^ + —(38462 - 95 - 27r^) 

+r'^{24 + 66r - 164r2 + 30r^ + 12r^) + 3r^{36 - 33r + 3r^ + 2r3)(ln(^-^) In(^) - Li2{-)) ] . 



(39) 



(40) 



r 



The behavior of these solutions at large r are 
115260 



equation of ^f^T'^ for A^^'' . The behavior is 



9o{r) 



67r2 - 139 ^ / 1 



1152r2 



, , 115262 -67r2 + 67 ^ / 1 

^^(") = 1152;:^ + ^1;:^ 



(41) 
(42) 



These solutions depend on the constants in the angular 
velocity of the electromagnetic fields 60,62, which must 
be determined by another boundary condition that we 
do not know. The total energy flux of electromagnetic 
fleld is 



,(4) 



(44) 



No matter how we choose 60,62 which must be deter- 
mined by a boundary condition, we cannot make A^^"^ 
0(l/r) as the second solution at large r. This result 
implies that this solution cannot work well at large r be- 
cause perturbation method breaks down there. However, 
the evaluation of the total energy flux is valid because, 
as mentioned above, E is independent of r whether we 
use perturbation method or not. 



E 



-2tt 



24 



7r(56 - 377^) 



2^2 "K^^-^" J 4^2 

1080 



(43) 



This is the value we want to evaluate. Fortunately, this 
result does not depend on the undetermined constants 
60,62. Probably, it is by chance. On the other hand, at 
large r, we can solve the r component of the fourth order 



III. COMPARISON WITH NUMERICAL 
CALCULATION 

In order to sec how well the forth-order term works 
to describe the total energy flux for large Kerr param- 
eter, we have performed numerical simulations of the 
monopole solution and compared the results with the an- 
alytical, second-order and forth-order solutions. 
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We have developed a two-dimensional General Rela- 
tivistic Ma gne to- Hydro Dynamics (GRMHD) code fol- 
lowing [ll| [l2l|. We have adopted a conservative, shock- 
capturing scheme with Harten, Lax, and van Leer (HLL) 
flux term ^13i] with flux-interpolated constrained trans- 
port technique We use a third-order TVD Runge- 
Kutta method for evolution in time, while Monotonized 
central slope-limited linear interpolation method is used 
for second-order accuracy in space [l5|. 2D scheme (2- 
dimensional Newton-Raphson method) is adopted for 
transforming between conserved variables and primitive 
variables [12l |. We used a simple gamma-law equation 
of state with 7 = 4/3. Modified Kerr-Schild coordinate 
is adopted with mass of the black hole (M) fixed where 
the Kerr-Schild radius r is replaced by the logarithmic 
radial coordinate xi = Inr. In the following, we use 
G = M = c=l unit. 

The computational domain is axisymmetric, with a 
grid that extends from nn — 0.98r-(- to rout — 230 
and from 6 = to 9 ~ tt where r_|_ is the outer event 
horizon. The numerical resolution is 300 x 300. As 
an initial condition, we put the 0th order terms of the 
monopole solution around the black hole [l6|. That is, 
sRf^ = -n*Ft"' = (0,aC sin 61/7=^, 0,0) in the Kerr- 
Schild coordinate where Ui,, *Ff^'^,a,g are the normal 
observer's four-velocity, the dual field tensor. Lapse- 
function, and determinant of the Kerr-Schild metric. 
The numerical constant C is set to be unity in this 
study. The plasma velocity relative to the FIDO is set 
to zero initially, and its pressure and density are set to 
the same value of P = p = Sft^/lOO so that force free 
approximation is a good one. Also, to keep the mag- 
netization reasonably low, when the critical condition 
0.013fJ2 > F^p -I- (7F2 - (7 -!)){/ is satisfied, density 
and internal energy U is increased by the same factor so 
that the critical condition holds [l^. Here F is the bulk 
Lorentz factor of the fiuid measured in the Kerr-Schild 
coordinate. We have performed numerical simulations 
with the Kerr parameters 0, 0.01, 0.05, 0.1, 0.15, 0.2, 
0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, 0.99, and 0.995 until 
time T = 200. 

In Fig. [2 we plot the total energy fiux at the final stage 
for small Kerr parameter (0 < a < 0.2) by rectangular 
points. The total energy fiux was evaluated at r = 20, al- 
though wc found that the total energy fiux is insensitive 
to the radius where it is evaluated. This means that the 
conservation of the total energy fiux has been confirmed 
numerically. Dashed line is just the interpolation of the 
calculated values. For comparison, the second-order ana- 
lytical solution is shown by the dotted line and the forth- 
order analytical solution is shown by the solid line. From 
this comparison, we can see that all of them coincide with 
each other that means the Blandford-Znajek solution is 
a really good approximation for small Kerr parameters. 

The situation becomes different drastically for large 
Kerr parameter. In Fig. [21 we plot the same values with 
Fig. [2 but with wide range of the Kerr parameter (0 < 
a < 1). We can see clearly the difference among three 
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FIG. 1: Comparison of the derived, conserved, total energy 
flux. Dashed line with rectangular points is numerical result 
for small Kerr parameter (0 < a < 0.2), dotted line shows the 
second-order analytical solution, and solid line represents the 
forth-order analytical solution. 



cases. However, we would like to emphasize the following 
point: we could confirm that the forth order analytical 
solution is a better approximation rather than the second 
order analytical solution. In fact, the deviation between 
the numerical result and forth order analytical solution is 
less than that between the numerical result and second 
order analytical solution. Of course, we could also see 
that much higher-order analytical solution is required to 
reproduce the numerical result, which is out of the scope 
in this study. 



IV. DISCUSSION AND CONCLUSION 

A. discussion 

We have tried to solve the fourth order terms for 
the Blandford-Znajek monopole solution to evaluate ex- 
tracted energy fiux more accurately. Since the equation 
for B'^ is a first order differential equation, we could solve 
by imposing a boundary condition at the horizon. How- 
ever, Atj, obeys the second order differential equation, 
thus we need to put two boundary conditions. Bland- 
ford and Znajek (1977) imposed the boundary conditions 
at the horizon and infinity. The boundary condition at 
infinity was chosen to connect the Michel's solution. In 
this study, we tried to impose the same boundary con- 
ditions for the fourth order solution. However, we found 
that the fourth order solution can not be connected to 
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FIG. 2: Same with Fig.l, but for wide range of the Kerr 
parameter (0 < a < 1) 

Michel's solution with any boundary condition at infin- 
ity. After all, we could not find the fourth order solution 
of A,p assuming the force free condition. Furthermore, 
we found the perturbation method for Kerr parameter a 
breaks down at large r because of the behavior A^^"* ^ . 
Fortunately, the fourth order term of the total energy flux 
extracted from a rotating black hole, which we want to 
evaluate most, does not depend on the solution we could 
not solve. 

The boundary condition which we could not specify 
was considered in [13, [11] ■ In general, the equation for 
the stream function which is called Grad-Shafranov 
equation has some singular surfaces [l^ [13] ■ In the force 
free limit, this equation is singular at horizon and at 
light cylinder. Thus, in order to determine the solution 
uniquely, we have to impose proper boundary conditions 
at horizon and light cylinder. However, boundary condi- 
tions to be imposed is not understood well |17i] . 

We have compared the numerical result with our fourth 
order solution. The numerical result we use is the result 
in GRMHD simulation. But for the evaluation of the 
total energy flux, GRMHD result is almost same as the 
result in GRFFE (General Relativistic Force Free Elec- 
trodynamics) simulation [2l|, so we would get the same 
conclusion when we compare the GRFFE result with our 
fourth order solution. However, for the details of the 
electromagnetic fields (e.g. shape of the field lines), the 
numerical result of GRMHD would be different from the 
one in GRFFE. As our future work, we want to examine 



the difference among our analytical result, the numerical 
result of GRMHD, and the one of GRFFE. 

In our calculation, we assume an infinitely thin disk, 
and cannot sclf-consistently determine a constant C in 
Eq. (lisp which determines a magnitude of the magnetic 
field. However, in a realistic accretion disk, C may de- 
pend on a 22] and total energy flux may be determined 
self-consistently. Thus, when Blandford-Znajek mecha- 
nism is considered with a realistic accretion disk, our 
result Eq. (|15)) . which also contains the constant C, has 
to be treated carefully. 



B. conclusion 

Our aim of this paper is to evaluate the total en- 
ergy flux extracted from a rapidly rotating black hole 
by Blandford-Znajek mechanism more accurately. As a 
result, although we could not the obtain the all pertur- 
bation solution up to the fourth order, we could evalu- 
ate the total energy flux extracted from a rotating black 
hole without any ambiguity irrespective of the unknown 
boundary condition at infinity. Also, we found the per- 
turbation method for the fourth order terms breaks down 
at large r by solving the equation for A^^'^ under the force 
free condition. This would be because in the monopole 
solutions, the force free condition can not be compatible 
with a rotating black hole at the fourth order of Kerr 
parameter. 

From the comparison between the numerical solution 
that is valid for < a < 1 and the fourth order solution, 
we find that the fourth order solution reproduces the nu- 
merical result better than the second order solution. At 
the same time, since the fourth order solution does not 
match well with numerical result at large a, we conclude 
that more higher order terms are required to reproduce 
the numerical result. 
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